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velocimetry , was further modified to provide additional laser 
and pressure probe access. The purpose of the research was to 
devise laser anemometry measurement techniques in a confined 
annulus and improve the ability to compare numerical 
predictions with experimental results. Flowfield computations 
were completed using a viscous flow solver, with the numerical 
exit plane coincident with experimental measurement location. 
A data reduction program was written to transfer non- 
dimensional numerical output to the experimental coordinate 
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and flow angle. Endwall two-dimensional laser measurements 
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I. 



INTRODUCTION 



The non-intrusive measurement of secondary flows in 
annular turbines can provide insight into the design and 
performance of turbomachinery. Turbine stage efficiency can 
be improved with more precise understanding of the flow 
mechanics and the ability to numerically predict the response 
to blade redesign without manufacturing and testing. 

Goldman and Seasholtz [Refs. 1-3] have conducted 
continuing research in the field of laser-Doppler velocimetry 
(LDV) measurement and numerical prediction in an annular 
turbine cascade. Their efforts have concentrated on obtaining 
comparisons ahead, in, and downstream of the blade passage. 
In order to obtain their results, an extensive optical system 
with correction for curved surfaces has been developed. To 
obtain two-dimensional velocity, they have measured one 
component ard then rotated the optics to resolve the second 
component . 

This report presents the results of laser-Doppler 
velocimetry techniques and three-hole pressure probe 
measurements in the wake of a high turning turbine nozzle. 
Only radial fiber-optic probe traverses were performed. 
Circumferential surveys were obtained by manually rotating the 
annular cascade blade row upstream of the measuring location. 
Pressure measurements were compared to numerical results 
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obtained using a three-dimensional viscous computer program. 
Laser-veloci metry techniques using two-component fiber-optic 
equipment ware developed for the small annulus cascade 
including two-dimensional measurements of the endwall flow 
through a small access hole in the outer casing. 
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II. 



EXPERIMENTAL SETUP 



A. TEST FACILITY 

Airflow for the annular turbine cascade experiment was 
provided by a VA-312 Allis-Chalmers 12 stage axial-flow 
compressor housed in the Naval Postgraduate School 
Turbopropulsion Laboratory. The compressor was operated at 
12,000 rpm at various discharge pressures. This air supply 
was also used to generate up to Mach 0.795 through a 107.95 mm 
(4.25 in) diameter nozzle for calibration of a three-hole 
pressure probe. 

Upstream air pressure was metered through the use of two 
dump valves and one cutoff valve prior to a plenum chamber. 
Air from the plenum was directed to the test section through 
flow straighteners in a 254 millimeter (10 in) steel pipe as 
shown in Figure 1. Flow stagnation pressure was measured 
upstream at two locations. One combination probe was 
connected to a mercury manometer board located near the dump 
valve controls. This pressure-setting probe also provided a 
digital readout of flow stagnation temperature. The second 
probe, and the static ports in the air supply pipe, were 
connected to a Scanivalve for automated pressure acquisition. 

LDV seed material, atomized-glycerin particles, was 
introduced through an eight millimeter (5/16 in) diameter 
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copper pipe as shown in Figure 1. Seed atomization was 
accomplished using laboratory compressed air with a commercial 
TSI, Inc., Six-Jet Atomizer. Laser traverse table buffet from 
the exiting flow was eliminated with a prefabricated aluminum 
cone (Figure 1) . Flow characteristics were examined with the 
cone on and off, using pressure measurements and dye injection 
from the twc chord inner hub static ports. 




Figure 1. Supply Pipe Schematic 

B. ANNULAR TURBINE CASCADE 
1. GENERAL DESCRIPTION 

The Annular Turbine Cascade (ATC) was developed to 
provide a small-scale device for flow experiment and 
computation on a transonic turbine stator. Thomas [Ref. 4] 
developed the basic apparatus. Figure 2 shows a cutaway view 
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of the construction of the ATC with a radial view of the blade 
passage just upstream of the test section. Figure 3 shows the 
blade geometry on a section perpendicular to the radial plane 
and Figure 4 shows an end view with the wake positioning 
capability. The inner hub radius was 98.81 mm (3.89 in) with 
constant blade geometry to the outer radius of 116.46 mm 
(4.585 in). The stator row had 31 blades with a midspan 
spacing of 21.77 mm (0.857 in). Reference 4 contains the 
original set of drawings used to manufacture the ATC. 
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Figure 4. Wake Positioning System 
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2. PRESSURE MEASUREMENT MODIFICATIONS 



The original ATC pressure measurement access was 
located two axial chords downstream of the blade trailing 
edges. Figure 5 details the old and new probe holes. The 
outer casing was modified to accommodate the half-chord access 
hole and is shown in conjunction with the adjustable holder. 
A new set of static pressure ports was added to the inner hub 
as shown. The ports were set 90 degrees apart and in axial 
alignment with the original holes. Tubing was connected to 
provide an average pressure tap for both hub axial locations. 
Appendix A contains a complete set of the drawings used for 
machining. 
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3 . LDV MODIFICATIONS 



LDV access was changed to minimize light reflection 
and improve off-axis options. The desire to ultimately obtain 
a third component necessitated an offset angle for a second 
fiber-optic probe. The design allowed for simple removal and 
replacement of the window during operation for cleaning, and 
radial access as shown in Figure 6. The position of the wake 
relative to the probe volume was controlled externally by 
rotating the blade row, inlet bellmouth and center body ahead 
of the measurement volume. An alignment blank was used to 
orient the laser traverse axis to the ATC. The 1.0922 mm 
(0.043 in) diameter hole in the blank provided an endwall flow 
measurement capability as shown in Figure 7. 




Figure 6. LDV Access Modifications 
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C. PRESSURE PROBE 

Total pressure, Mach number and flow angle measurements 
were made with a three hole yaw probe [Ref 4]. All pressure 
taps were connected to a Scanivalve with a ±344 kPa (±50 psi) 
differential transducer. Scanivalve control and data 
recording were accomplished using an HP-9000 computer system. 
Appendix B contains calibration curves and a listing of the 
Hewlett-Packard code used to acquire and output data. 
Pressure data from the ATC were acquired with the same 
hardware and software used in the probe calibration. 

Previously, yaw alignment was accomplished using a "U"- 
tube water manometer. Difficulty encountered with hysteresis 
and the response of the water to yaw misalignment was 
corrected wi’zh the use of a differential pressure transducer. 
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D. 



LASER-DOPPLER VELOCIMETER 



The LDV system is shown schematically in Figure 8. The 
Lexel Model 95 four-Watt Argon-Ion LASER was connected to a 
TSI, Inc., Model 9201 ColorBurst Multicolor Beam Separator. 
The beam separator split the incoming light into two beams, 
with one beam passing through a Bragg cell for 4 0 Mhz 
frequency shifting. The two beams were further split into 
three polarized pairs: green (514.5 nm) , blue (488 nm) , and 
violet (476.5). The TSI, Inc., Model 9230 ColorLink 
Multicolor Receiver provided the frequency-shift signal to the 
ColorBurst for beam shifting. Each beam was directed through 
individual couplers for transmission to the laser probe by 
fiber-optic cable. The fiber-optic probes contained receiving 
optics to direct the return signals to the ColorLink. 




Figure 8 . LDV Wiring Schematic 
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The ColorLink provided photomultiplier and frequency 
shifting functions. Control of shifting and photomultiplier 
voltage was accomplished remotely by computer. After 
conditioned by the ColorLink, the signals were sent to 
individual TSI, Inc., Model 1990C signal processors where 
valid Doppler signals were counted with reference to desired 
characteristics set by the user. An MI-990 interface 
monitored coincidence criteria (2-D measurements) and 
transferred valid counter information to computer memory. 

The fiber-optic probes were mounted to a Model 41583 
LINTECH Traverse Table. An Applied Motion Products System 
1618 Traverse Controller was used to control traverse movement 
manually and by computer. [Ref. 5] All remote control of LDV 
system equipment and data processing was completed through the 
menu driven environment of TSI^s FIND version 4.04. 

E . EXPERIMENTAL PROCEDURE 
1. PROBE CALIBRATION 

Probe calibration and investigation of data 
acquisition options were completed together. The calibration 
was accomplished using a free-jet with a 107.9 mm (4.25 in) 
nozzle connected to the air supply (Figure 1) used for the 
annular cascade. Figure 9 shows the testing apparatus and 
pressures used to calibrate the cobra probe. 

The probe was aligned in yaw with a ±68.9 kPa (±10 
psi) Statham, SCO 277-10, differential pressure transducer. 
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Yaw-alignment response was an improvement over the previous 
manometer arrangement. The transducer bridge output was 
periodically zeroed by removing the probe from the flow. The 
yaw pressure measurement was tapped from one side of the 
transducer as shown in Figure 9 . The Scanivalve system 
recorded the five pressures shown in Figure 9. Pre-operation 
tubing-leak checks were completed by trapping pressure to each 
Scanivalve port in turn. 




Figu:;e 9. Yaw Probe Calibration Equipment 



The probe was calibrated at Mach numbers from 0.20 to 
0.79 in steps of 0.1. Pitch sensitivity was investigated at 
twelve positions from minus ten degrees to plus eight degrees 
(Figure 9) . Appendix B contains a description of the 
calibration variables and shows the sensitivity of the probe 
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to pitch angle. Sixth-order polynomials were fit to the 
calibration data. The zero-degree probe calibration was 
successfully validated (Appendix B) under different 
atmospheric conditions. 

2. PRESSURE PROBE MEASUREMENTS 

Wakes were surveyed by repositioning the blade passage 
upstream relative to the probe. The inner hub rotated with 
the blade ring as the wake was positioned for LDV or pressure 
measurement. The bolt shown in Figure 4 was used to position 
and hold the inner ring and stator. A total of 18 degrees 
circumferential variation, in steps of one degree, was 
achieved. 

Cobra-probe measurements were completed as described 
in the calibration section. Upstream static and stagnation, 
and inner-hub static pressures were used in conjunction with 
the center-port and yaw-port pressures of the probe. The 
probe provided total pressure and flow angle. Mach number was 
obtained from (3 (Appendix B) . The desired hub-static to 
upstream-stagnation pressure ratio was set using the 
Scanivalve system and manual control of the air supply dump 
valves. The probe was positioned radially and aligned with 
the flow before pressure data were acquired. Upstream 
stagnation temperature and wake flow angle were recorded 
manually for each measurement location. Flow angle was read 
from a vernier scale to an accuracy of plus or minus 0.2 
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degrees at span locations greater than 25 percent. Span 
position was established with reference to 100 percent at the 
outer wall. 

3 . LDV ALIGNMENT 

LDV alignment for endwall flow measurement was 
attempted as shown in Figure 10. The laser mounting bracket 
micrometer was set to ten millimeters for visual alignment 
around the access hole as shown. Obtaining a repeatable depth 
reference was not possible. This alignment technique did 
provide micrometer scale coverage for data recording and 
ensured the beams passed cleanly through the hole. The length 
of the probe volume precluded the use of the outer surface of 
the plate as a reference point. 
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Alignment for automatic traverse of the wake passage 
through the window was accomplished using the alignment blank 
(Figure A8) . The multi-component beams were positioned in 
each of the six alignment holes using manual control of the 
traverse table. The coordinates were read from the digital 
position indicator on the traverse controller for each hole. 
Axis conversion from the traverse to the ATC was completed 
using the position data and the dimensions of the alignment 
blank (Appendix A) . Traverse table movement could then be 
programmed using the automatic traverse control menu in TSI's 
FIND software. 

4. ENDWALL MEASUREMENTS 

Endwall flow data were acquired without the use of 
automatic traverse control. The laser beams were aligned with 
the half-axial-chord access hole as described above, and then 
traversed manually using the probe-bracket micrometer. 
Adjustments to seeding position, counter gain and 
photomultiplier voltages were continually updated to maximize 
data rate and minimize the effect of noise. Seed-wand 
position was particularly important. Micrometer position 
information was married to the appropriate data files in post 
processing of FIND's output ASCII .PRN files using Borland 
Corporation's spreadsheet program Quattro Pro version 5.0. 
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III. COMPUTATIONAL FLUID DYNAMICS 
A. GRID GENERATION 

Grid generation was completed using Turbomachinery C GRID, 
TCGRID, [Ref . 6]. The code was developed to build 3-D "C” or 
"H" type grids for turbomachinery in the standard PL0T3D 
format (NASA Ames Research Center graphics software) . 
Appendix C contains the namelist input file used to generate 
the stator g.iometry and a discussion of the output. The inlet 
bellmouth and center body as shown in Figure 4 were not 
modeled for this study. 

Careful manipulation of the namelist variables was 
required to obtain a smooth solution. The grid geometry was 
developed to improve resolution of the endwall boundary layer. 
The ability to generate a fine grid was hampered severely by 
the length of wake desired in the model. Reduction of the 
wake-exit plane from one axial chord downstream of the blade 
trailing edges [Ref. 4] to one-half axial chord resulted in a 
dramatic increase in grid resolution. The final grid 
dimension was 150x31x65, up from 121x31x21. The blade surface 
received an additional 29 points and the radial resolution was 
improved from 21 points to 65. The smaller computational wake 
length also reduced the interval between grid lines from the 
trailing edge of the blade. 
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B . COMPUTATIONAL SCHEME 



Flow solutions were completed using Rotor Viscous Code 3-D 
(RVC3D) . RVC3D [Ref. 7] was designed to solve three- 
dimensional viscous flows in turbomachinery. The code solves 
thin-layer Navier-Stokes equations in Cartesian coordinates 
using an explicit finite difference technique. Turbulence 
effects were modeled using the Cebeci-Smith option. The 
restart option of the code was utilized throughout the 
investigation to monitor convergence of the solution. 
Appendix D contains the namelist input file used to obtain the 
final solution and a description of steps required to run the 
code on the Naval Postgraduate School (NPS) Cray computer. 
Output solution files (fort. 3) were examined using the 
graphics software PL0T3D. Comparison with experiment required 
data reduction of the output file. Appendix E contains a 
synopsis of the output file structure with data reduction 
steps required for comparisons and a listing of the FORTRAN 
program (plane. f) used to accomplish the manipulation. 

The final solution selected for comparison was chosen for 
the minimum residual and the smoother transition of a solution 
across the wake cut. An area of unsteady flow became 
prominent near the hub and on further analysis revealed an 
area of predicted reverse flow. Figure 11 contains a 
comparison between total pressure ratio results at the 50 
percent span location for 1200 and 1490 iteration solutions. 
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One blade wak^e is shown with the wake cut in the grid clearly 
evident in the plot. The closely-spaced points in the wake 
cut of the TCGRID geometry has introduced a mathematical 
boundary. The 1200 iteration solution exhibited the beginning 
of the reverse flow area near the hub, but contained serious 
discontinuities because of the solution scheme. The 1490 
iteration total-pressure-ratio plot demonstrates the smoothing 
that took place as the solution progressed. 



CFD Solution Comparison 

Total Pressure Ratio (50 % Span) 




Figure 11. CFD Solution Comparison, 50% Span 
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IV. RESULTS AND DISCUSSION 
A. PRESSURE PROBE 

Annular Turbine Cascade pressure measurements were 
obtained at an average Reynolds number of 650000 based on 
axial chord. The pressure ratio, inlet total to exit static 
at the hub, was set at 0.68 to enable a proper comparison to 
be made with the CFD solution. Pressure and flow angle 
measurements were repeatable, with Scanivalve pressure- 
acquisition standard deviations of less than 0.20 percent. 
However, data obtained inside 25 percent span were at low 
dynamic pressure and too difficult to resolve with the 
pressure probe, and are therefore not reported. Flow 
unsteadiness increased as the probe was immersed further in 
from the enc'wall. The direction of approach to zero yaw in 
the less-than-25-percent span region affected the flow angle 
measurement by as much as 0.5 degrees. 

Pressure and flow-angle data were taken at 18 wake 
positions for 25, 50, 75 and 90 percent span locations on the 
half-axial-chord wake measurement plane. Pressure data 
repeatability was demonstrated along the 50 percent span 
station for all wake positions, and radially at the two-chord 
zero-degree wake position. Figure 12 shows the total- 
pressure-ratio data obtained for the 50 percent span location. 
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The data show that there was a negligible influence of the 
conical diffuser. Figure 13 shows the repeatability obtained 
in radial surveys for a constant wake position at the two- 
chord measurement location. The repeatability shown in 
Figures 12 and 13 clearly validated the ability to reset the 
desired operating pressure ratio for the experiment. Flow 
angle measurements were repeated at the 90 percent span 
location, and radially at the two-chord wake measurement 
location as shown in Figures 14 and 15. In all cases, the 
repeatability of the probe pressure measurements was 
excellent, with the largest discrepancy less than 0.5 percent 
in total pressure. 
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Figure .12. Pt/PO Repeatability, 50% Span 
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Figure 13. Pt/PO Repeatability, Radial 
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Flow Angle Repeatability 

Half chord (90% Span) 




Wake Position (degrees) 



Figure 14. Angle Repeatability, 90% Span 



Flow Angle Repeatability 

2 Chord (Diffuser Ofi) 




Figure 15. Angle Repeatability, Radial 
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B . LASER MEASUREMENTS 



Two-dimensional LDV measurements were obtained through a 
1.0922 mm (0.043 in) hole to a depth of approximately three 
millimeters. Initial measurements were taken with one 
component (green) of the two-component probe. Figures 16 and 
17 show the velocity profile and turbulence intensity results 
from the first attempt at endwall measurements. The 'depth' 
shown in the graphs is an estimated distance from the endwall 
based on the rapid loss of data rate observed at the depth 
shown as 0.1 mm from the endwall. The data map a boundary 
layer profile to approximately 0.6 millimeters. The high 
pressure ratio (0.98) provided near perfect seed conditions, 
and the casing vibrations were low. The data shown for V-z do 
not necessarily coincide with the exact locations of the data 
for V-0 because the fiber-optic probe was rotated 90 degrees 
and alignment could not be maintained. 

The turbulence intensity results are with respect to the 
maximum velocity at the furthest measurement location into the 
annulus. The significantly higher turbulence intensity in the 
V-0 direction is probably due to the steep gradient of the V-0 
mean-flow profile in the radial direction. 
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LDV Endwall Flow Velocity 

Pressure Ratio - 0.98 
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Figure 16. Single Component LDV Endwall Velocity 



LDV Endwall Turbulence Intensity 
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Figure 17. Single Component LDV Turbulence 
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Two component LDV measurements were obtained at 
approximately 15 percent higher flow velocity. Figure 18 
depicts a test section 'sweep' of near-endwall velocity 
measurements at a pressure ratio of 0.974. Depth reference 
(in mms with respect to an unknown datum on a vernier scale) 
is shown. The exact distance from the wall was not 
obtainable. Flow periodicity is suggested and decreasing 
velocities toward the outer casing show the endwall boundary 
layer. The )>oundary layer is distorted due to corner vortices 
and secondary flows that originate in the blade row. Flow 
angle and turbulence intensity plots are included in Figures 
19 through 21. Once again the turbulence intensity in the V-0 
direction is higher than that in the V-z direction. Figures 
18 through 21 represent automatic editing of the histograms by 
the FIND software (user selection) . Data outside two standard 
deviations were discarded. 

Increased flow velocity introduced higher vibration levels 
and temperature increases as the test progressed. Hole 
alignment was periodically verified and, occasionally, small 
adjustments (less than 0.01 mm) were required. 
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Depth Reference (mm) 
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Figure 18, Two Component LDV Endwall Velocity 



LDV Endwall Flow Angle 

Depth Reference (mm) 




Wike Position (degrees) 



Figure 19. Two Component LDV Flow Angle 
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Vtheta Turbulence Intensity 

Referenced to Maximum Velocity 
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Figure 20. Two Component LDV Turbulence, V-theta 
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Figure 21. Two Component LDV Turbulence, V-z 
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C. NUMERICAL RESULTS 

The CFD solution was run with the same pressure ratio 
(0.68) used in the experiment. Graphical output of the 
numerical results was useful in obtaining a qualitative 
understanding of the flow characteristics. Figure 22 depicts 
the predicted pressure ratio for the blade surface at the hub, 
mid-span and tip. The suction surface curves suggest that the 
position of the blade passage throat was at 0.8 of axial 
chord. Figure 23 shows, with colored contours, pressure 
coefficient information for the blade surface pressure 
distribution. The Mach number contours shown in Figure 24 
provided an initial indication of where flow characteristics 
were well established. Boundary layers appear as expected, 
thin near the endwall and thick near the hub. Wake curvature 
in the spanwise direction is also evident in the Mach contour 
plot. Difficulty in obtaining repeatable flow angles in the 
experiment below 25 percent span may have been due to the 
unsteady flow separation which is suggested by the 
computationa.. results for the region near the hub wall. 
However, the reverse flow predicted by the code in this region 
could not be verified in the experiment. 
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Figure 22. RVC3D Blade Pressure Ratio Plot 
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Figure 23. PL0T3D Blade Pressure Coefficient Plot 
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Figure 24. PL0T3D Exit Plane Mach Contour Plot 
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D. COMPARIS-)N 



Three comparisons were made between numerical and 
experimental results. To show the flow periodicity , the CFD 
solutions were repeated over two blade passages as shown by 
the solid line in the following Figures 25 through 36. 
Comparisons between experiment and numerical (Appendix E) 
results are presented for total pressure ratio, Mach number 
and flow angle. Figures 25 through 28 show the total pressure 
ratio comparison. It is noted that uncertainty in the probe 
measurement of total pressure is small. The probe calibration 
uncertainty enters only when comparing Mach number. Wake- 
profile phase coincidence was established using the 50 percent 
span Mach number data. The phase was kept constant for the 
circumferential comparisons at all other spanwise locations. 
Mach number (Appendix B) comparisons are contained in Figures 
29 through '2, and Figures 3 3 through 3 6 contain the flow 
angle data. 

CFD total pressure ratio predictions were within 10 
percent of the experimental values at 25 percent span. The 
difference decreased to below one percent at 90 percent span. 
Mach number comparisons with measurements were within 15 
percent at 25 percent span and decreased to three percent at 
90 percent span. Flow angle comparisons were within five 
degrees at 90 percent span with improvement to within two 
degrees at 25 percent span. The reason for closer agreement 
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in total pressure ratio and Mach number toward the tip and 
improving agreement in flow angle toward the hub is 
unexplained. 
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Figure 25. Total Pressure Ratio Comparison, 25? 
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Figure 26. Total Pressure Ratio Comparison, 50% 
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Total Pressure Ratio 

Comparison (75 % Sp>an) 
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Figure 27. Total Pressure Ratio Comparison, 75% 
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Figure 28. Total Pressure Ratio Comparison, 90% 
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Mach Number 

Comparison (25 % Span) 
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Figure 29. Mach Number Comparison, 25% 
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Figure 30. Mach Number Comparison, 50% 
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Mach Number 

Comparison (75 % Span) 
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Figure 31. Mach Number Comparison, 75% 



Mach Number 

Comparison (90 % Span) 




OD • t«p I 



Figure 32. Mach Number Comparison, 90% 
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Wake Flow Angle 

Comparison (25 % Span) 
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Figure 33. Flow Angle Comparison, 25% 
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Figure 34. Flow Angle Comparison, 50% 



38 



Wake Flow Angle 

Comparison (75 % Span) 
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Figure 35. Flow Angle Comparison, 75% 
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Figure 36. Flow Angle Comparison, 90% 
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CONCLUSIONS AND RECOMMENDATIONS 



Experimental measurements and computational output 
locations were made to coincide. Repositioning the 
measurement volume to one-half axial chord downstream from the 
blade trailing edges improved the computational grid 
resolution. Flow solutions with the revised grid more 
realistically resolved endwall boundary layer characteristics 
so that future endwall LDV measurements will have a 
corresponding computed solution with equivalent resolution. 

Pressure probe surveys will continue to be important in 
comparing experiment with CFD prediction. The cobra probe 
calibration obtained was reliable and repeatable. An 
extensive set of baseline data has been accumulated for future 
use. The Scanivalve data acquisition system demonstrated the 
ability to obtain pressure readings with a standard deviation 
of less than 0.2 percent. The HP-9000 output was manually 
entered into a spreadsheet program for post processing. The 
time required to manually enter the output was overshadowed by 
the spreadsheet data manipulation flexibility. 

The most definitive comparison between CFD and probe 
measurements was in the total pressure ratio. Probe 
calibration error was not present in this comparison. The 
total press ire ratio disagreed by ten percent at the 25 
percent spar location and less than one percent at the 90 
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percent location. In all cases, the computer solution 
predicted a lower ratio. The difference may be attributed to 
computational modeling parameters or bias introduced in the 
experimental pressure ratio setting. Upstream total pressure 
was measured at only one station and this value may not have 
been a proper reflection of the average total pressure. Mach 
number comparison characteristics followed the trend of the 
total pressujre ratio as expected. Flow angle comparison with 
CFD was best at 25 percent span with a maximum difference of 
5° at the 9C percent span location. Of all the measurement 
locations, the 90 percent span location was the best defined 
and flow angle measurements made there were repeatable within 
0.2 degrees. 

Endwall LDV through a small (1.0922 mm) casing access hole 
has been demonstrated. The ability to resolve endwall 
boundary layers and eventually tip flow losses appears to be 
promising. Precise LDV probe-volume positioning techniques 
were not obtained. The maximum pressure ratio attained for 
endwall LDV was 0.974. CFD and pressure probe data were 
accumulated at an operating pressure ratio of 0.68. Low data 
rate, vibration of the endwall hole, and thermal expansion 
problems will have to be dealt with as the operating pressure 
ratio is decreased to 0.68. 

Future efforts with the computational fluid dynamics 
should include modeling the inlet bellmouth and center body. 
The parameters selected for use in RVC3D should also be 
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explored to obtain a qualitative feel for cause and effect on 
the solution with particular attention to the unsteady region 
of the exit flow. 

Pressure data acquisition can be continued to obtain a 
more complete map of the wake region. Pressure 
instrumentation of a blade and passage would provide an 
additional source of CFD code validation data. Data reduction 
and presentation would be improved with a change in operating 
system. The current Scanivalve system with the HP-9000 does 
not allow file transfer to other computers for post 
processing. HP-9000 limitations restrict output flexibility. 

The annular turbine cascade now has off-angle LDV options 
that have not been explored. Three-dimensional radial LDV 
surveys of the wake passage should be attempted for comparison 
with pressure measurements already obtained. Initial 
indications show that the LDV probe-volume energy may not be 
able to provide return signals of sufficient intensity as the 
flow velocity increases. The addition of corrective optics to 
the fiber-optic probe may be required. 

The endwall technique developed in this study did not 
provide pressure equalization across the access hole. The 
pressure across the access hole should eventually be equalized 
to eliminate flow through the hole and therefore any influence 
the opening might have on the flow structure [Ref. 9]. The 
hole acted as a natural optical noise filter and further 
investigation of hole diameter influence on LDV performance is 
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warranted. In order to truly understand the hole influence 
and measurement reliability, the glycerin seed particle 
distribution and a study of their momentum characteristics 
should be undertaken. With these refinements, access to the 
flow through a hole in the casing is feasible, without 
changing the endwall flow characteristics. 
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APPENDIX A. DRAWINGS 



Figure Al. 



O 

T5 



<D 




Adjustable Pressure Probe Holder 
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Figure A2 . ATC Mounting Flange Modification 
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Figure A3. Outer Ring Modifications 
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Figure A4 
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Outer Ring Modification, Section A 



47 




48 






Figure A6 
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Outer Ring Modification, Section C 
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APPENDIX B. COBRA PROBE CALIBRATION 



The following synopsis was extracted from Appendix A of 
Reference 8. Mach number or dimensionless velocity can be 
related to a non-dimensional pressure coefficient, f3 , as 
follows; 



Total velocity: 






V,= 



Dimensionless velocity: 



x= 






Dimensionless velocity and Mach number are related as follows; 



c.T, 






M^fRT 

so that 



= . 2 



i+X_1m2 

2 



or 



, 2 - 2^2 






(y-l) il-X^) 
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A pressure coefficient for the probe aligned in yaw can be 
written as 



where p^ and p^ are probe sensor measurements. Using p^ and P23 
respectively, it can be shown that 



in the range of velocities which are of interest. 

Figure Cl contains the calibration data without curve 
fitting and Figure C 2 shows the validation data plotted 
against the sixth order and linear cuirve fits at zero pitch 
angle. The sixth order curve fit equation is 




-P23 =f[M, Pitch) 
Pi 



Z=(a<9+aip+a2p2+a3p3+a4p^+a5p5+a5p®)xl00000 



with the coefficients shown in Figure Bl. 
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Figure Bl. Cobra Probe Calibration Coefficients 
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Figure B2. Cobra Probe Calibration Data Plot 
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Figure B3 . Probe Validation 
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Data acquisition was accomplished using 



an HP-9000 



computer system with the following program: 



10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

271 

272 
280 
290 
300 

310 

311 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 



!FIIE NAME : C0BRA_ATC 

!DISK LABEL: "/AUSTIN" 

! UT JEFF AUSTIN 

!IAST MODIFIED 02/03/94 BY LCDR JOE SPITZ 

I 

! THIS PROGRAM RECORDS AND REDUCES MEASURED PRESSURES 
! FRCM A Scanivalve OCM®CrED TO A THREE HOIE PRESSURE 
! PROBE. Scanivalve PORT OCMTOCITCWS ARE AS FOLLOWS: 

I 

! PROGRAM VARIABLES AND PORT OCM®Cn(»LS. . . (#3 not used) 

! V = DESIRED S/V 

! A = PRESENT S/V PORT 

] C = CHANNEL 

! VO = DVM READING 

! S$ = SCANNER LISTEN CODE 

! S = SCANNER # 

! Scanivalve #5 SET TO ADVANCE CMLY 

! Pamb = Ambient I^ressure (psi) 

! Patm = Atmosi±ieric pressure. . .Pt>rt # 1 

! Peal = (Calibration pressure. . .Pt>rt # 2 

! Pt.ot = Reference total pressure. . .Ptort # 4 

! Phub = Inner hub static pressure. . .Pt>rt # 5 

! Pstat = Upstream static pressure. . .Port #6 

! PI = (Center hole on prc±>e. . .Port # 7 

! P23 = Yaw port pressure. . .Pt>rt # 8 

! beta = (P1-P23)/P1 

! prat = Phub/Ptot (pressure ratio) 

! mach = mach(Pl,Phub,(Sainma) 

I 

OFTTCN BASE 1 
(Camma=1.4 
CLEAR SCREEN 
PRINTER IS (CRT 

DISP "PIEASE WATT WHILE RESETTING ScanivalveS" 

PRINT 

V=5 

A=1 

QDSUB Read 

J *****★*>?**************************★★*************★**★*★*■*★*★**★** 
PRINT "DATA ACCQUISITTON FOR THREE HOLE YAW PROBE" 

DIM Y(20) 

DIM Pat3Ti(10,10) , Peal ( 10, 10) ,Ptot( 10,10) , Phub (10,10) , Pstat (10, 10) 
DIM PI (10, 10), P23(10,10) 

I 

INPUT "ENTER MONTH, DAY, YEAR (IE. 02,05,94)" ,Y(3) , Y(5) , Y(7) 
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431 INFUT "ENTER RUN #: ", Y(9) 

432 FRINTER IS 711 

440 PRim? USING "K,DD, DD, DD","DftTE OF RUN: ",Y(3) , Y(5) , Y(7) 

441 PRINT "DATA RUN ",Y(9) 

450 PRINT 

451 PRINTER IS CRT 

460 !***************************************************************** 
550 PRINT 

560 PRINT "CmSTANT FLOW. WAKE: +9 TO -9 DEGREES. SPAN: 5 IDCATIONS 

570 PRINT 

571 PRINT " Zero DVM on Scanivalve #5, port #1." 

572 PRINT 

580 PRINT " Set 20 inHg for Scanivalve calibration." 

581 PRINT 

582 PRINT" Adjiast span for DVM = 0.009823 volts on port #2." 

600 PRINT 

610 PRINT" ENSURE DVM ADJUSTED AND DESIRED PRESSURE RATIO SET. . ." 

630 FOR WAKE = 1 TO 10 

640 INPUT "ENTER AMBIENT PRESSURE IN iriHg",Pinhg 

641 Paitib=Pinhg* ( . 4911541) 

650 PRINT 

660 PRINT "CHECK DESIRED SPAN AND FIOW AITO^MENT SET. . ." 

661 PRINT 

662 PRINT "*****IRESS OCOTTNUE WHEN READY TO TAKE DATA*****" 

670 PAUSE 

671 CLEAR SCREEN 

680 PRINT 

700 FOR SEAN = 1 TO 10 

710 V=5 

720 A=1 

730 G0SU3 Read 

740 ********tv***************************************************** 

750 1 

840 1 ************************************************************* 

850 IREOORD IHE PRESSURE DATA FOR THIS WAKE AND SPAN 
860 1 
870 ! 

880 
890 
900 
910 ! 

912 

913 

914 

915 

916 

919 

920 

921 ! 

922 1 



OUTPUT 722;"FIR1TIM0Z1 " ! SETS UP HP3456A DVM 

V=5 ! Scanivalve #5 DESIGNATED 

S=1 ! HP3495A SCANNER #1 DESIGNATED 

PRINTER IS 711 
PRINT 

PRINT "Wake Data Groi^) ", WAKE 
PRINT "Span Location Counter ", SEAN 
PRINT USING "/,5X,K,DDD,/,2X,K,15X,K,17X,K";" 
Scanivalve #",V, "Port", "Volts", "Psi" 

PRINT 

PRINTER IS CRT 

READ Scanivalve PORTS 
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923 ! 

930 

940 

950 

960 

970 

980 

990 

1000 

1002 

1010 

1020 

1021 

1022 

1023 

1024 

1025 

1026 
1027 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 

1040 

1041 

1042 

1043 

1044 

1045 

1046 

1047 

1050 

1051 

1052 

1053 

1054 

1055 

1056 



FOR A=1 TO 8 
GOSUB Read 
WATT 2.0 

OtmUT 701 ;"C" ! CLEARS SCANNER #1 

OOTIVr 701 USING ''DCO'';V+9 J SCANNER SET TO READ Scanivalve 
PORT 
TRIGGER 722 
ENTER 722 ;V0 

PRINTER IS 711 
! Pressure in psi 
IF 1^1 THEN 

Patm (WAKE , SPAN) =VO*1000 

FRENT USING ”2X,MX),7X,7DD.DDDDDC®,8X,7E».IXXX)"; 

P, VO, Patm (WAKE, SPAN) 

FRENTER IS CRT 

FRENT USING ''2X,IXO,7X,7rX).DDCmX),8X,7CO.E»DD''; 

P, VO , Patm (WAKE , SPAN) 

IRENTER IS 711 
END IF 

J 

IF A=2 THEN 

Peal (WAKE , SPAN) =VO*1000 

FRINT USING "2X,E»D,7X,7r».DD[XXXXD,8X,7rX).DEXXl"; 

P, VO , Peal (WAKE , SPAN) 

FRENTER IS CRT 

FRINT USING ''2X,rX»,7X,7CO.DDKXXX),8X,7Ct).[XX»"; 

P, VO, Peal (WAKE, SPAN) 

FRENTER IS 711 



IF A=4 THEN 

Ptot (WAKE , SEAN) = (V0*1000) -Patm (WAKE, SPAN) +Pamb 
FRINT USING "2X,DDD,7X,7CO.rXXXm3,8X,7Ct).DDtXi"; 

P, VO , PtCJt (WAKE , SPAN) 

FRENTER IS CRT 

FRINT USING ''2X,DDD,7X,7CO.EO[XX»D,8X,7IX).DDDD'' ; 

P, VO , Ptot (WAKE , SEAN) 

FRENTER IS 711 
END IF 

IF A=5 THEN 

Riub (WAKE , SPAN) = ( V0*1000) -latm (WAKE , SPAN) +Pamb 
FRINT USING ''2X,IXX),7X,7[X).EOXm3,8X,7DD.DDE»"; 

P, VO , FHub (WAKE , SPAN) 

FRENTER IS CFO? 

FRINT USING "2X,IX»,7X,7C».IXXXXXDD,8X,7CO.DDDD" ; 

P, VO , Fhub (WAKE, SEAN) 

PRINTER IS 711 
END IF 
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1057 

1060 

1061 

1062 

1063 

1064 

1065 

1066 

1067 

1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 

1080 

1081 

1082 

1083 

1084 

1085 

1086 

1087 

1088 

1089 

1090 

1091 

1092 

1097 

1098 

1099 

1100 

1101 

1102 

1103 

1104 

1105 

1106 

1107 

1108 

1109 

1110 

1111 



I 



IF A=6 THEN 

Pstat (WAKE , SPAN) = (VO* 1000) -Patm (WAKE , SPAN) +Pairib 
FRD7T USING "2X,DCO,7X,7[XD.rXDDIXt)D,8X,7DD.IXXX3''; 

P , VO , Pstat (WAKE , SPAN) 

BRINTER IS CRT 

FKENT USING "2X,DC»,7X,7IX).IXt)CmD,8X,7DD.EODD''; 

P, VO , Pstat (WAKE , SPAN) 

HRINTER IS 711 
END IF 

I 

IF A=7 THEN 

PI (WAKE , SPAN) = (V0*1000) -Patra (WAKE , SPAN) +Pamb 
PKCNT USING ”2X,DC»,7X,7DD.DDCOE»D,8X,7DD.CODD"; 

P , VO , PI (WAKE , SPAN) 

PKCNTER IS CRT 

FKENT USING ''2X,DCO,7X,7DD.DCCDCOD,8X,7IX).C®C»"; 

P,V0,P1 (WAKE, SPAN) 

FKENTER IS 711 
END IF 

I 

IF A=8 THEN 

P23 (WAKE, SPAN) = ( VO* 1000 )-Patati (WAKE, SPAN) +PaiTib 
FRENT USING ”2X,IXX),7X,7DD.DDDDCDD,8X,7DD.EXXO" ; 

P,VO, P23 (WAKE, SPAN) 

FKENTER IS CRT 

FKENT USING ''2X,DDD,7X,7DD.EXX)DC»D,8X,7DD.!X03"; 

P , VO , P2 3 (WAKE , SPAN) 

FKQTTER IS 711 

I 

prat=Fhub (WAKE , SPAN) /Ptot (WAKE , SPAN) 
beta= (PI (WAKE , SPAN) -P23 (WAKE, SPAN) ) /PI (WAKE, SPAN) 
inach=SQRT (5* (PI (WAKE , SPAN) /Ftiub (WAKE , SPAN) ) -^ ( 1/3 . 5) -5) 
END IF 

I 

NEXT A 

FKENT "prat = ", prat 

FKENT "beta = ", beta 

FKENT '•mach # = ", mach 

FKENTER IS CRT 

FKENT 

FKENT "Set Span and eilign prcbe. . .CONTTNEJE WHEN READY. . . " 

FKENT 

PAUSE 

CLEAR SCRE04 
NEXT SPAN 

FKENT "Verify Wake & i^xdate Panto. . .OCS'lTIt^UE WHEN READY. . . " 

FKENT 

PAUSE 

CLEAR SCREEN 
NEXT WAKE 
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1112 

1113 

1114 

1115 

1116 

1117 

1118 

1119 

1120 

1121 

1122 

1123 

1124 

1125 

1126 

1127 

1128 

1129 

1130 

1131 



I *****************Erro OF MAIN PROGRAM*********************** 

I 

! SUBROOTINE TO POSITICX^ AND READ S/V PORTS 
! Scanivalve #1 HAS HOP HOME FEATURE 

Read:OUTHJr 707 USING '#,K";V 
P0=SPOLL(707) 

L=BINAND(P0,15) 

Tl=RCrCATE(P0,4) 

M1=BINAND(T1,7) 

P=10*M1+L 
CIEAR 707 

IF P=A THEN Finish 
OUTFUr 701; "C" 
oumir 701 USING ''IXO";V-l 
OUIRTT 701; "C" 

WATT .1 
GCOO Read 
Finish: RETURN 
END 
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APPENDIX C. GRID GENERATION 



Grid generation was completed on IRIS workstations. 
TCGRID, a FORTRAN program, required the adjustment of 
parameter statements to dimension arrays for the model build. 
Parameter adjustment required patience. The author found the 
manual helpful in most cases, however, the most useful tool 
turned out to be a running summary of the namelist input for 
troubleshooting. The most difficulty was observed in 
obtaining a smooth solution on the suction side of the blade 
geometry. The following TCGRID namelist file was utilized to 
obtain the grid shown in Figures Cl and C2 . Reference 4 
contains the original list of blade coordinates and is not 
included here. 

&naml im=150 jro=31 km=65 itl=20 icap=18 k2d=3 merid=0 Send 
&nam2 nle=16 nte=10 dsle=.018 dste=.003 d^ub=.0004 dstip=.0004 
dswte=.001 dswej?=.060 dsthr=l. dsminF=.0004 dsmax=.025 
dsi'a=.45 rcom=.098 Send 

&nam3 itenrt=100 idbg=0 0 0 0 0 0 0 0 aahto^.S codd=.45 Send 
&nam4 zbc=-1.5 -1.5 0.5 -1.5 -1.5 .5 

rbc^ 3.895 3.895 3.895 4.585 4.585 4.585 Send 
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Figure C2 . Grid Geometry, 3-D Perspective 
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APPENDIX D 



RVC3D CODE INPUT 



RVC3D runs were initially completed on IRIS workstations 

at NFS. Solutions were limited to less than 500 iterations to 

ensure residuals were decreasing. File transfer of the fort.l 

(grid) and the rvc3d.in (namelist) files to the Cray were 

completed only after initial testing was completed. The 

following namelist contains the commands used for the final 

set of iterations to the minimum residual attainable. It is 

a restart (iresti=l) from 1200 iterations with termination at 

1490 iterations (imax=290) . 

'TRANSONIC TORBINE Annular Cascade - (Restart 12-1490 ) ' 

&nll im=150 jm=31 km=65 itl=20 iil=67 &end 

&nl2 cl=5.5 aviscl=0.0 avisc2=0.0 avisc4=0.50 ivdt=l 

nstg=4 itmax=290 irs=l ^i=0.50 ^j=0.60 ^k=0.60 Send 
&nl3 ibcin=3 ibcex=3 isyrat=0 ires=10 icmt=50 
iresti=l iresto=l ibcpw=0 iqirFO Send 
&nl4 emxx=0.13 eratyO.O emrz=0.0 expt=0.0 prat=0.6800 
ga=1.4 om=0. 000000 igecm=l alex=-67.0 &end 
&nl5 ilt=3 tw=1.00 renr^.651e6 pmr^.7 prtr=.9 

vispwr=. 666666 srtip=0.0 crautin=14. jedgef=31 kedge=20 
iltirr=2 dbll^.0024 dblt=0.0089 Send 
&nl6 iol=l io2=165 oarK). ixjb=0 njo=l nko=3 jo=l 
ko=5 11 16 &end 

Runs on the NFS cray were accomplished with the command 
>qsub myjob.nqs 

where the file myjob.nqs was as follows: 
myjob.nqs 

#QSUB -q reg -IT 43200 -IM 9Mw 
cd cfd 

rvc3d <rvc3d.in> rvc3d.out 
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The RVC3D executable, grid (fort.l), restart (fort. 2), 
and the input namelist (rvcSd.in) were held in the 
subdirectory cfd. The final solution required 13.2 hours of 
Cray computer time. A plot of the residuals is shown in 
Figure Dl. 



RVC3D Residuals 

Normalized to Maximum 




Figure Dl. RVC3D Residuals 
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APPENDIX E 



RVC3D DATA REDUCTION 



Exit plane data were reduced using the FORTRAN program 
plane. f. RVC3D output the following non-dimensional values 
for each grid point in the file fort. 3: 



<7i 



<?2 






<?4 



<?5 



Po 

=_£iL 

Po<^o 

= _£21 
Po^o 

pw 

Po'^0 



=_EL 



Cn^ 2 1 Cn 






w 



2 M 



/J 



Computational output was manipulated in plane. f for 
comparison with experimental results. The following table 
provides a coordinate system comparison. 



Plot3d 


Annular 


Cascade | 


, Cartesian 


Cartesian 


cylindrical p 


x(i) - u 


z(i) - q2 


z(i) - q2 ! 


y(i) - V 


y(i) - q3 


V«(i) 


z(i) - w 


x(i) - q4 
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Transformation to cylindrical coordinates was conducted 
using the grid location information contained in the RVC3D 
grid input file fort.l. 

0(i)=sin~^f 1 z{i) =y/yTTy^^i^F(iY^ 

All angles (6) were referenced to the y axis shown in the 
blade geometry plot of Figure 3 . Flow angle computations were 
completed using the non-dimensional q values described above 
as follows: 

Flow Angle T (dimensional) ; 

Ve = - (u) sin6+ ( v) COS0 V^= (u) cos0+ ( v) sin0 V^=w 

Y^tan-if-^ 



Flow Angle 



V'e =-^^ = - ig^) sin0+ (gj) cos0 
Po^o 



F, =-^=g2 



^nd n C' 






Vq ^ 



^ad/ 



Direct comparison was possible with pressure ratio 
measurements. The ratio between local, and upstream reference 
total pressures was computed as follows: 



-£.=y(y-i) 

Po 



<? 5 " 



2 Po 




+ 









M 



/J 
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p{u^ + v^ + w^] 



Po Po 2po 



^=Y(Y-i)<?5 + 




Mach number was obtained using the isentropic 
relationship between Mach number and temperature ratio. 



The program plane. f is included here for reference. It 
is an updated version of the program written by Thomas [Ref. 
4]. The nev; version includes a calculation of flow angles 
after conversion to cylindrical coordinates for proper 
comparison with experimental results. 



c********************************************************************* 

c 

c plane. f reads rvc3d files & writes ascii files for plotting 
c exit plane mach number and flow angles 
c (ccsrparison with ejqperimental results) 

c unit 1 = irpat xyz file 

c unit 3 = iiput q file 

c unit 4 = output residual history 

c unit 7 = output (k. #1) counter, theta, pt/pO, Mach, Angle 
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c unit 8 = output (k #2) cxjunter, theta, pt/pO, Mach, Angle 

c unit 9 = output (k #3) counter, theta, pt/pO, Mach, Angle 

c unit 10 = output (k #4) counter, theta, pt/pO, Mach, Angle 

c unit 11 = output (k #5) counter, theta, pt/pO, Mach, Angle 

c 

Q** *********************** ******** ************************** ********** 

parameter (ni=150 , nj =3 1 , nk=65 ) 
integer urn 

real x(ni,nj,nk) ,y(ni,nj,nk) ,z(ni,nj,nk) , gama 

qq(5,ni,nj ,nk,resd(5000,5) ,ql,q2,q3,q4,q5,m(5) 
dimension kk(5) ,ang(5) ,theta(5) ,tt0(5) ,vtnd(5) ,vmd(5) , 

1 vznd(5) ,pt^(5) ,deg(5) 
c 

c k-valuss are hard-wired belcw (10%, 25%, 50%, 75%, 90%) 
c (use span.f to chtain %-k equivalency) 

c 

data kV21,27,33,39,45/ 

Q** ******************************************************* ************ 

c read q.rid coordinates 

c********************************************************************* 
read(l,*) im, jm,km 

read(l,*) (((x(i,j,k) ,i=l,im) ,j=l,jm) ,k^l,km) , 

1 (((y(i,j,k) ,i=l,im),j=l,jm),k^l,km) , 

2 (((z(i,j,k) ,i=l,im) ,j=l,jm) ,k=l,km) 

Q** ********************************************************* ********** 

c read restart file 

c********************************************************************* 

read(3,*) imax, jmax,kmax 
read(3,*) fsmach,aljha,re,titne 
c 

icheck=iabs ( im-imax) +iabs ( jm- jmax) +iabs (km-kmax) 
i f ( icheck . ne . 0 ) then 
write(6,610) im, jm,km,imax, jmax,kmax 
stop 
endif 
c 

read(3,*) ( ( ( (qg(l,i, j,k) ,i=l,im) ,j=l,jm) ,k?=l,km) ,1=1,5) 
c 

c additional residual data 

read (3,*) itl,iil,fhdeg,gama,om,nres,dum,dum,dum,dum 
read(3,*) ( (resd(nr,l) ,nr=l,nres) ,1=1,5) 

c********************************************************************* 

c Pressure side of exit wake cut. . . 

c***********^-********************************************************* 

k=kk(3) 

i=l 

do 10 j=jmax, 1,-1 
um=unH-l 
do 15 1=1,5 
k=kk(l) 
c 



70 



c 

c 



Rename desired q values for use 



ql=qq(l,i,j,k) 

q2=gq(2,i,j,k) 

q3=qq(3,i,j,k) 

q4=qq(4,i,j,k) 

q5=qq(5,i,j,k) 

c 

c Ccxipute pt/pO ratio 
c 

ptpO (1) =gama* (gama-1) *q5+(gama-gama**2/2) /ql* (cp**2+q3**24q4**2) 
c 

c Conpute T/Tt ratio and mach number 
c 

ttO(l)=gama* (gama-l)/ql*(q5-0.5* (q2**2+q3**2+q4**2)/ql) 
m ( 1 ) = ( ( 1/ttO ( 1 ) -1) *2 . 0/ (gama-1 ) ) * *0 . 5 
c 

c Ccsrpute Iheta and Non-D velocities theta, r and z 
c 

theta(l)=asin(y(i,j,k)/(z(i,j,k)**2+y(i,j,k) **2)**0.5) 
vtnd(l)=-q4*sin (theta (1) )+q3*cos (theta (1) ) 
vmd(l)=q4*cos (theta (1) )+q3*sin (theta (1) ) 
vznd(l)=q2 

deg (l)=theta (1) *57 . 296 
c 

c Conpute flow angle (degrees) using vtnd and vznd 
c 

ang(l)=atan2(vtnd(l) , vznd (1) ) *57.296 
15 continue 

write(7,300)um,deg(l) ,ptp0(l) ,m(l) ,ang(l) 
write (8, 300) urn, deg (2) ,ptp0(2) ,m(2) ,ang(2) 
write(9,300)um,deg(3) ,ptp0(3) ,m(3) ,ang(3) 
write(10,300)um,deg(4) ,ptp0(4) ,m(4) ,ang(4) 
write (11, 300) urn, deg (5) ,ptp0(5) ,m(5) ,ang(5) 

10 continue 
c 

c Repeat above on other side of wake cut at imax 
c j=l not used (duplicate location across wake) 

c jmax not repeated (same location as j=31 at i=l) 

c 

jmaxiTf=j max-1 
i=im 

do 30 j=2,jmaxm,l 
umFum+1 
do 25 1=1,5 
k=kk(l) 

ql=qq(l,i,j,k) 
q2=qq(2,i,j,k) 
q3=qq(3,i,j,k) 
q4=qq(4,i, j,k) 
q5=qq(5,i,j,k) 
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ptpO(l)=gama* (gaina-1) *q5f (gam-gairB**2/2)/ql* (q2**2+q3**2+q4**2) 
tt0(l)=gama*(gaina-l)/ql*(q5-0.5*(q2**2+q3**2+q4**2)/ql) 
m(l) = ( (l/ttO(l)-l) *2/(gaina-l) ) **0.5 

theta(l)=asin(y(i,j,k)/(z(i,j,k)**2+y(i,j,k)**2)**0.5) 
vtnd(l)=-q4*sin (theta (1) )+q3*oos (theta (1) ) 
vmd(l)=q4*cx5s (theta (1) )+q3*sin(theta(l) ) 
vznd(l)=q2 

deg (1) =theta (1) *57 . 296 
ang(l)=atan2 (vtnd(l) ,vznd(l) ) *57.296 
25 cxjntinue 

write (7, 300) urn, deg (1) ,ptp0(l) ,ro(l) ,ang(l) 
writers, 300) urn, deg (2) ,ptp0(2) ,m(2) ,ang(2) 
write(9,300)um,deg(3) ,ptp0(3) ,m(3) ,ang(3) 
write(10,300)um,deg(4) ,ptp0(4) ,ro(4) ,ang(4) 
write (11, 300) urn, deg (5) ,ptp0(5) ,m(5) ,ang(5) 

30 continue 

c residual history output to unit 4 

c********************************************************************* 
write(4,310) 1, (resd(l,l) ,1=1,5) 
do 40 j=2,nres 
it=10*(j-l) 

40 write (4, 310) it, (resd(j ,1) ,1=1,5) 

300 format(i5,5f8.3) 

310 format(i5,5(l>c,el0.3)) 

610 format(' ***** warning *****',/, 

1 ' iin, jm, km, read frcm iiput',3i5,' do not match',/, 

2 ' im, jm, km, read from restart file',3i5) 
stop 

end 
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